Correlated Dirac Fermions on the Honeycomb Lattice 
studied within Cluster Dynamical Mean Field Theory 
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The role of non-local Coulomb correlations in the honeycomb lattice is investigated within 
cluster dynamical mean field theory combined with finite-temperature exact diagonalization. 
The paramagnetic semi-metal to insulator transition is found to be in excellent agreement with 
finite-size determinantal Quantum Monte Carlo simulations and with cluster dynamical mean field 
calculations based on the continuous-time Quantum Monte Carlo approach. As expected, the 
critical Coulomb energy is much lower than within a local or single-site formulation. Short-range 
correlations are shown to give rise to a pseudogap and concomitant non-Fermi-liquid behavior 
within a narrow range below the Mott transition. 



I. INTRODUCTION 

The recent discovery of graphenei has greatly stimu- 
lated the study of the electronic properties of the honey- 
comb latticed In view of the vanishing density of states 
at the Fermi level, an issue of particular interest is the 
influence of electron-electron interactions. Gonzalez et 
al£ performed renormalization group calculations and 
showed that the suppression of screening of the long- 
range Coulomb interaction gives rise to deviations from 
conventional Fermi-liquid behavior. Lattice field the- 
ory simulations 4 indicated a Coulomb driven second- 
order semi-metal to insulator transition. Meng et al& 
performed extensive variational Quantum Monte Carlo 
(QMC) simulations for the Hubbard model with varying 
cluster sizes and identified a spin-liquid phase between 
the semi-metallic state characterized by massless Dirac 
fermions and an antifcrromagnetically ordered Mott in- 
sulator. The onset of the long-range antiferromagnetic 
order was found to be consistent with previous QMC cal- 
culations for finite-size clusters The Mott transition 
of the honeycomb lattice was also investigated^ within 
single-site dynamical mean field theory! (DMFT). How- 
ever, because of the small number of nearest neighbors, 
the neglect of spatial correlations in this system is ques- 
tionable and gives rise to a significant overestimate of 
the range of the semi-metallic behavior up to large val- 
ues of the onsite Coulomb interaction. To account for 
non-local correlations in the honeycomb lattice, Wu et 
alJ^ recently applied a cluster extension^ of DMFT 
(CDMFT) by using continuous-time QMCi^ as impu- 
rity solver. The transition between the semi-metallic 
and insulating phases was found to occur at a consid- 
erably smaller critical Coulomb energy than within the 
single-site DMFT and to be in good agreement with the 
variational QMC results by Meng et al^ 

In the present work we use finite-temperature ex- 
act diagonalization 14 (ED) in combination with clus- 
ter DMFT to investigate the two-dimensional Hubbard 
model on the honeycomb lattice for unit cells consisting 
of six sites. The focus is on the dynamical properties 
of the non-local self-energy which have not been studied 



before. Moreover, in view of the large size of this unit 
cell and the approximate nature of quantum impurity 
solvers, CDMFT results obtained within complementary 
schemes are clearly desirable. An important advantage of 
ED is the accessibility of large Coulomb energies and low 
temperatures, and the absence of sign problems. Also, 
in contrast to finite-size variational QMC, ED is appli- 
cable away from half-filling. On the other hand, due to 
the exponential growth of the Hilbert space, the number 
of levels representing the bath surrounding the cluster 
is severely limited. Here, we use 12 levels in total, i.e., 
six impurity levels and six bath levels. Since these bath 
states are coupled indirectly via the onsite Coulomb re- 
pulsion within the six atom cluster, the spacing between 
excitation energies is very small. Finite-size errors are 
thereby greatly reduced, even at low temperatures, so 
that self-energies and spectral functions can be evaluated 
reliably at rather low real frequencies. 

The ED/CDMFT results discussed below reveal a con- 
tinuous Mott transition in excellent correspondence with 
the variational QMC simulations by Meng et aL— and 
with the QMC/CDMFT calculations by Wu et alH The 
critical onsite Coulomb energy is considerably smaller 
than the one found in single-site DMFT calculations^ 
Furthermore, short-range correlations included within 
CDMFT are shown to give rise to metallic and insulat- 
ing contributions to the self-energy at the Dirac points 
in the Brillouin Zone, where the former dominate at low 
Coulomb interactions, and the latter increase essentially 
quadratically with the nearest-neighbor non-local self- 
energy component. These terms lead to the excitation 
gap above the Mott transition. Below the transition they 
yield a narrow pseudogap. Thus, short-range correlations 
induced via onsite Coulomb repulsion give rise to devia- 
tions from Fermi-liquid behavior in some range below the 
critical interaction strength. Also, the effective mass en- 
hancement does not diverge at the Mott transition, but 
increases to a finite value. The opening of the pseudogap 
below the transition, and the variation of the effective 
mass with Coulomb energy, are qualitatively similar to 
analogous results obtained within cluster DMFT calcu- 
lations for the square lattice^ - — 
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The outline of this paper is as follows: In the next 
section we briefly outline the application of ED / CDMFT 
to the Hubbard model for the honeycomb lattice. Section 
III provides the discussion of the results. The summary is 
presented in section IV. The focus in this work is on the 
paramagnetic semi-metal to insulator transition. Spin 
liquid and antiferromagnetic phases will be addressed in 
a subsequent publication. 



II. CLUSTER ED/DMFT FOR HONEYCOMB 
LATTICE 



<55 



0.4 



0.3 



0.2 



0.1 



~1 1 1 T" 




In this section we discuss the combination of cluster 
DMFT with finite-temperature ED for the purpose of 
evaluating the effect of non-local Coulomb interactions 
on the honeycomb lattice. The Hubbard Hamiltonian is 
given by 

h = -t ]T ( c +c ia + h.c.) + uJ2 mtmi, (!) 

where the sum in the first term includes only near- 
est neighbors, t is the hopping matrix element, and U 
is the onsite Coulomb repulsion. The band dispersion 
for the honeycomb lattice may be written as e(k) = 
±i|l + e jfclV ^ + e l ( fexV ^+' c !' 3 )/ 2 |. In the following we define 
t = 1 as energy unit. 

Let us divide the two-dimensional lattice into clusters 
consisting of six sites. Within the unit cell, the positions 
are specified as ai = (0,0), a 2 = (1,0), a 3 = (v / 3/2,3/2), 
a 4 = (V3,l), a 5 = (V3,0), and a 6 = (^3/2, -1/2). 
The nearest neighbor spacing is taken to be a = 1. 
The supercell lattice vectors are then given by A 1 / 2 — 
(3\/3/2,±3/2). Within CDMFT— the interacting lat- 
tice Green's function in the cluster site basis is given by 

Gij(iu> n ) = ^ i ioJ n +fJ>- t(k) - E(iuj n )]T. , (2) 

k 

where uj n = (2n + l)7rT are Matsubara frequencies and 
H is the chemical potential. The k sum extends over the 
reduced Brillouin Zone, t(k) denotes the hopping matrix 
for the superlattice, and Ey (jw n ) represents the cluster 
self-energy matrix in the site representation. The diag- 
onal elements of the symmetric matrix G^ are identical 
and there are three independent off-diagonal elements: 
G12 = Gi6, G13 = G15 and G14. Because of these sym- 
metry properties, it is convenient to go over to a diagonal 
"molecular orbital basis" , in which the elements G m (iLv n ) 
are given by 

G1/2 = (G11 + 2G13) ± (G14 + 2G12) 
G3/4 = (G11 — G13) ± (G14 — G12) (3) 
= G 5 / 6 . 

The self-energy matrix satisfies the same symmetry 
properties as G and can therefore be diagonalized in 



FIG. 1: (Color online) Total density of states p(^) (solid 
curve) of honeycomb lattice and cluster components p m (w) 
(dashed curves) within diagonal molecular orbital basis. For 
clarity, these components are divided by n c = 6. Orbitals 3 
and 4 are doubly degenerate, u — defines the Fermi energy 
for half-filling. 

the same fashion. These elements will be denoted as 
£ m (zu; n ). Below we focus on the special case of half- 
filling. Since the density of states is then particle-hole 
symmetric with respect to uj = 0, Gu(iui n ) is purely imag- 
inary. The same applies to G13, whereas G12 and G14 are 
real, corresponding to odd density of states components. 
Thus, the diagonal molecular orbital components of G 
satisfy G 2 = -G\ and G 4 = —G%. Figure [1] illustrates 
the uncorrelated density of states components in the di- 
agonal basis, where p m {to) = -jImG m (w) for S = 0. 
The even and odd onsite and intersite components may 
be obtained by inverting Eq. (j4|. 

A central feature of DMFT is that, to avoid double- 
counting of Coulomb interactions in the quantum impu- 
rity calculation, the self-energy must be removed from 
the small cluster in which correlations are treated explic- 
itly. This removal yields the impurity Green's function 

G (lLU n ) = [GilLUn)- 1 + E(iLJ n )]-\ (4) 

For the purpose of perfoming the ED calculation we now 
project the diagonal components of Go(iw„) onto those 
of a larger cluster consisting of six impurity levels and 
six bath levels, i.e., n s = 12 is the total number of levels. 
Thus, 

G , m (iu) n ) ~ G C Q m {iuj n ) 

= \iw n + fi - e m - . Vmk — I ,(5) 

where e m denotes impurity levels, bath levels, and V m k 
hybridization matrix elements. The incorporation of the 
impurity level e m in the fitting procedure yields a more 
accurate representation of Go^ m (iu} n ) than by projecting 
only onto bath levels. 

Assuming independent baths for the diagonal cluster 
molecular orbitals, each component Go >m (iw„) is fitted 
using three parameters: one impurity level e m , a bath 
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level ek and a hopping integral V m k- To evaluate the 
finite-temperature interacting Green's function of the 
cluster it is useful to transform the impurity orbitals back 
to the site representation in which the Coulomb interac- 
tion is diagonal. We denote this transformation by T, 
where the matrix elements are given by 
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(6) 



Thus, the diagonal 6x6 subblock of the cluster Hamilto- 
nian, hf, = (efe^fefe')) representing the bath levels remains 
unchanged, while the 6x6 impurity subblock becomes 
nondiagonal in the cluster site basis: 
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with 
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r' 



[(ej + e 2 ) + 2(e 3 + e 4 )]/6 
[(ex - e a ) - (e 3 - e 4 )]/6 
[(ei + ea) - (e 3 + e 4 )]/6 
- e 2 ) + 2(e 3 - e 4 )]/6. 



(8) 



We point out that the hopping element t of the origi- 
nal lattice Hamiltonian does not appear since it is ab- 
sorbed into t via the molecular orbital cluster levels 
e m which are adjusted to fit Go tm (icu n ). The procedure 
above therefore includes not only hopping between clus- 
ter and bath. It also introduces four new parameters 
within the six-site cluster: the onsite level e, and up to 
third-neighbor hopping parameters: r, r', and r". At 
half-filling, 62 = — ei and e 4 = —£3 for symmetry rea- 
sons, so that e = t' = 0. In this mixed site- molecular 
orbital basis, the hybridization matrix elements V m k be- 
tween cluster and bath molecular orbitals introduced in 
Eq. ([5]) are transformed to new hybridization matrix el- 
ements between cluster sites i and bath orbitals k. They 
are given by 



V- k = {TV) ik =Y,T im Vmk ■ 



(9) 



The single-particle part of the cluster Hamiltonian now 
reads 



(10) 



Adding the onsite Coulomb interactions to this Hamilto- 
nian, the non-diagonal interacting cluster Green's func- 
tion at finite T can be derived from the expressio n 18 ' 19 
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(11) 



where E v and \ v) denote the eigenvalues and eigenvectors 
of the Hamiltonian, j3 — l/T and Z = J2 V exp(— jSE^) is 
the partition function. Further details concerning the 
evaluation of the cluster Green's function can be found 
in Ref . [2C] where the analogous procedure is discussed for 
the square lattice. Since G^J satisfies the same symme- 
try properties as G^ , it is diagonal within the molecular 
orbital basis, with elements G^. The diagonal cluster 
self-energy components are then given by an expression 
analogous to Eq. (U): 



(12) 



The important assumption in DMFT is now that this 
impurity cluster self-energy is a physically reasonable 
representation of the lattice self-energy. Thus, 
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where, at real frequencies, T, m {oj) is continuous whereas 
S^(w) is discrete. 

In the next iteration step, these diagonal self-energy 
components are used as input in the lattice Green's func- 
tion Eq. @ , which in the molecular orbital basis may be 
written as 



G m {iuj n ) 



k 



[iu) n + fi — Tt(k)T 1 - E(iw„)] 



(14) 

where T is the transformation defined in Eq. ((6]). Note 
that Ti(k)T~ 1 is not diagonal at general k points. As a 
result, all molecular orbital components of S(iw„) con- 
tribute to all components G m (iuj n ). We also point out 
that, to get adequate resolution at low frequencies, be- 
cause of the vanishing density of states, a sufficiently 
large number of k values near the Dirac points must be 
included in the Brillouin Zone integration. 



III. RESULTS AND DISCUSSION 

Figure 2(a) shows the low-energy region of the inter- 
acting density of states for several Coulomb energies, at 
temperature T — 0.005. These distributions are derived 
from an extrapolation of the local lattice Green's function 
Gn(iuj n ) to real frequencies. To illustrate the stability 
of this extrapolation, at each value of U several curves 
are plotted for 50 to 200 Matsubara points, with an ad- 
ditional small energy broadening of the order of O.lw 2 . 
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FIG. 2: (Color online) (a) Low-energy region of density of 
states An(w) = — - Im Gn (ui) of honeycomb lattice for sev- 
eral Coulomb energies at T = 0.005. The noninteracting den- 
sity of states is indicated by the black dotted curve. Between 
50 and 200 Matsubara points are used to extrapolate the lat- 
tice Green's function to real frequencies, (b) Density of states 
over wider energy range for U = 5 and U = 7. 



(For > 1 the broadening is kept constant at 0.1.) At 
U = 3, a tiny gap or pseudogap is seen which is near 
the limit of what can be resolved within ED/DMFT. At 
U = 4, a full gap of width A w 0.25 has opened. Its width 
increases approximately to A « 0.6 when the Coulomb 
energy is increased to U = 5. This trend is consistent 
with the one found in Refs. [Mill. The variation of the 
gap over a wider range of U is indicated in Fig. 3(a). 
The spectral distributions in Fig. 2 show that the van 
Hove singularity at ui = ±1 is strongly broadened and its 
weight is shifted to lower energies. Above the transition, 
the Hubbard bands are difficult to resolve as long as U is 
less than the band width, but they become pronounced 
once U > W, as shown in Fig. 2(b) for U = 7. 

These results suggest that non-local correlations in the 
honeycomb lattice induce a paramagnetic semi-metal to 
insulator Mott transition in the range U = 3 ... 4. Be- 
cause of the continuous nature of the transition (see be- 
low), it is difficult to identify the precise value of the 
critical interaction. Nevertheless, our finding is con- 
sistent with the variational QMC simulations 5 ' and the 
QMC/DMFT calculations^ which yield U c « 3.6. It is 
also in qualitative agreement with earlier finite-size clus- 
ter QMC simulations which gave U c ~ 4.5 (Ref. @) and 
17 C «4...5 (Ref.0). On the other hand, all of these val- 



(c) 




0.2 



v> o 



13 



T = 0.005 



14 
12 



3.5 



4 
U 



4.5 




(d) 



FIG. 3: (Color online) (a) Excitation gap A, (b) average dou- 
ble occupancy d occ , (c) local and nonlocal spin correlations 
(SizSiz), and (d) molecular orbital occupancies n m as func- 
tions of Coulomb energy for T = 0.005. Orbitals 3 and 4 
are doubly degenerate. There is no indication of hysteresis 
behavior in the critical region U = 3 ... 4. 



ues are significantly lower than the ones obtained within 
single-site DMFT which yields (Xw 10 . . . 13 £2, More- 
over, in agreement with Refs. [a-lWllI we do not find any 
hysteresis behavior for increasing versus decreasing U, 
as shown in Figure 3(b) for the double occupancy, in- 
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dicating that the transition is continuous. In contrast, 
within local DMFT the transition was shown to be of 
first order £ Fig. 3(c) shows the onsite and intersite spin 
correlations, (S\ z Si Z ), for i = 1 . . .4. The onsite and sec- 
ond neighbor components are positive, while the first and 
third neighbor components are negative, underlining the 
antiferromagnetic nature of the spin correlations. 

One of the interesting effects of Coulomb interactions 
in multi-orbital systems is the possibility of correlation- 
induced charge transfer between orbitals. As shown in 
Figure 1, the six-site unit cell of the honeycomb lat- 
tice maybe viewed as consisting of six molecular orbitals 
which are split by an effective crystal field and there- 
fore have different orbital occupancies. Figure 3(d) shows 
the variation of these occupancies with Coulomb energy. 
Evidently, there is little orbital polarization, a result 
that was also observed in CDMFT calculations for the 
square and triangular lattices 3i Moreover, the double 
occupancy, the spin correlations, and the orbital occu- 
pancies reveal no clear sign of a Mott transition in the 
region where the spectral distribution exhibits the open- 
ing of a gap. 

To analyze the nature of the semi-metal to insulator 
transition, it is therefore necessary to examine the non- 
local contributions to the self-energy. Figure 4 shows the 
four independent components of the cluster self-energy 
S(iw n ) within the site basis, for Coulomb energies in the 
region of interest, U — 3 ... 5. For symmetry reasons, 
Eii(iw„) and T,i^(iuj n ) are purely imaginary. They be- 
have as ~ iu) n at low frequencies. In contrast, Y^i2{iu) n ) 
and £i4(«cj„) are real and approach a finite value in the 
limit uj n = 0. 

In a seminal paper long before the synthesis of 
graphene, Gonzalez et al& studied the influence of 
electron-electron interactions on the quasiparticle life- 
time in a single layer of graphite. Taking into account 
the long-range nature of the Coulomb interaction, their 
renormalization group calculations indicate that the sup- 
pression of electronic screening at low frequencies yields 
deviations from conventional Fermi-liquid behavior, with 
ImE(w) approximately linear in uj for uj — 0.4. . .3 eV 
(for t = 2.4 eV). 

To determine possible non-Fermi-liquid contribu- 
tions to the self-energy derived within the present 
ED/CDMFT approach, we have carefully searched 
for w„ln(u; n ) behavior in the imaginary components 
En(iw„) and Si3(jw„). Within the accuracy of our re- 
sults, these functions do not indicate any such deviations 
and seem to be well proportional to iuj n in the entire 
range U — ... 5. Also, they do not indicate a finite lim- 
iting value for uj n — > which would imply a finite lifetime 
for states near the Fermi energy. Thus, the non-Fermi- 
liquid properties obtained in Ref. [H seem to be associ- 
ated with the long-range part of the Coulomb repulsion 
which is absent in the Hubbard model for purely onsite 
interactions^ 

We note, however, that to understand the spectral 
features of the quasiparticle density of states, it is not 
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FIG. 4: (Color online) Self-energy components Eii(icj n ), i — 
1 ... 4, for honeycomb lattice in cluster site representation as 
functions of Matsubara frequency for Coulomb energies U = 
3 ... 5 in steps of 0.2; T = 0.005. 



sufficient to study the self-energy components shown in 
Fig. 4. In particular, these isolated components do not 
provide any evidence for a Mott transition in the region 
U = 3 ... 4, where the density of states shown in Fig. 2 in- 
dicates the opening of a gap. To illustrate the smoothness 
of the self-energy components in this range of Coulomb 
energies, we show in Fig. 5(a) the slopes of ImEu and 
ImEi3, and the values of Re £12 and Re £14 in the low- 
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FIG. 5: (Color online) (a) Low-frequency limits of self-energy 
components: slopes of ImSu(w„), Im£i3(icd n ), and values 
of Re E12 (ikin), Re£i4(iu; n ), as functions of U. (b) Amplitude 
b 2 /(l — a) of insulating contribution to T,(K, iu) n ), Eq. ((20}, 
as a function of U. T = 0.005 



frequency limit. Evidently, these individual components 
do not reveal the existence of the Mott transition seen 
in the density of states. This behavior differs qualita- 
tively from the Hubbard model for the square lattice at 
half-filling, where at the metal insulator transition the 
(-7T, 0) component of the self-energy changes from ~ iu n 
to ~ l/iui n at small uj n , and the real part of the (0,0) 
and (ir, tt) components exhibits a jumpi 16 i 

The origin of this apparent discrepancy is the fact that, 
as pointed out above, the local interacting density of 
states depends in a highly nonlinear manner on all non- 
local self-energy components £y. This is evident from 
the expression for the lattice Green's function, Eq. ©, 
where the hopping matrix t(k) and E(iw„) cannot be si- 
multaneously diagonalized, as indicated also in Eq. (Till) . 
To account for this admixture of intcrsite self-energy ele- 
ments, it is useful to examine the 6x6 cumulant matrix 

M(iu n ) = [iu n - £(kj„)] _1 . (15) 

Since M has the same symmetry properties as E, its non- 
local components are given by 

Mn = [(Mi+M 2 ) + 2(M 3 + M 4 )]/6 

Mia = [(Mi - M 2 ) - (M 3 - M 4 )]/6 

Mia = [(Mi+M 2 )-(M 3 + M 4 )]/6 

M14 = [(Mi - M 2 ) + 2(M 3 - M 4 )]/6 (16) 




FIG. 6: (Color online) (a) Imaginary part of self-energy at K, 
Eq. (|20[) , as a function of Matsubara frequency, for U = 1 ... 5 
at T — 0.005. (b) uj n lm'E(K,iuj rl ), demonstrating the range 
of the insulating part of the self-energy at Dirac points K. 
Solid red curves denote integer values of U, dashed blue curves 
denote intermediate values of. 



where the diagonal molecular orbital elements are 

M m (iuj n ) = [iuj n - E m (iu; n )] _1 . (17) 

The opening of the Mott gap takes place at the six 
K points of the Brillouin Zone. To analyze the behav- 
ior of the cumulant at these points, we make use of the 
periodization2£ 

1 6 

M(k, iu n ) =qJ2 e^-'^My (*u n ), (18) 

ij=l 

where aj are the positions within the six-site cluster. At 
K = 2tt(2/3v / 3, 0) and K' = 2tt(1/3v / 3, 1/3), this ex- 
pression simplifies to 

M(K,iu n ) = Mn(iuj n ) - M 13 (iuj n ) 

= [M 3 (iw n )+M 4 (iWn)]/2. (19) 

The self-energy at K is therefore given by 

E(K,iu> n ) = iio n - M~ l {K, iu> n ) 

sa iu n a+- — -, u„4 0, (20) 

iuj n {l - a) 

where a is the initial slope of Im [En(iui n ) — Ei 3 (iw n )] 
and b the low-frequency limit of Re [Ei 2 (io; rl ) — Ei 4 (io; n )]. 



7 



This self-energy is shown in Figure 6 for various Coulomb 
energies. The above expression indicates that H(K,iu) n ) 
is imaginary as expected for particle-hole symmetry at 
the Dirac points. It consists of metallic (~ iui n ) and in- 
sulating (~ l/iu; n ) contributions. The insulating term, 
which is responsible for the opening of the Mott gap, 
increases quadratically with b = Re [S12 — S14]. Thus, 
the semi-metal to insulator transition is driven primarily 
by the nearest-neighbor component of the non-local self- 
energy, with a minor additional contribution due to the 
third-neighbor self-energy, and a weak renormalization 
related to the initial slope of Im. [En(iui n ) — Ei3(iu; n )]. 
The variation of the amplitude b 2 /(l — a) of the insulat- 
ing term with Coulomb energy is depicted in Fig. 5(b). 
The comparison with Fig. 3(a) demonstrates that the ex- 
citation gap A roughly tracks the amplitude of this term. 

According to the results shown in Figure 4, a w 
-0.2 . . . -0.6 and b « -0.05 ... - 0.6 in the range U = 
3 ... 5. Thus at U = 3 the amplitude of the insulat- 
ing term is about 10 2 times smaller than at U = 5. 
Nevertheless, this small contribution is responsible for 
the pseudogap below the Mott transition, indicating the 
breakdown of Fermi-liquid behavior in the metallic phase. 
For U < 2 we find |6| < uiq, so that the pseudogap 
can no longer be resolved within the accuracy of ED. 
A similar pseudogap induced by short-range correlations 
at half-filling was observed below the Mott transition in 
the Hubbard model for the square lattice*^— Neglect- 
ing the small insulating term of iui n ) sufficiently far 
below the transition, the effective mass enhancement of 
the quasiparticle bands near the Dirac points is given by 
m*/m = 1 - a = 1.0 ... 1.2 for U = ... 3, i.e., it does 
not diverge at the Mott transition, in contrast to results 
derived within single-site DMFT. This finding is also con- 
sistent with the behavior seen on the square latticed 



IV. SUMMARY 

The influence of onsite Coulomb interactions on the 
electronic properties of the honeycomb lattice has been 
investigated within cluster dynamical mean field theory 
combined with exact diagonalization. The interacting 
density of states exhibits the opening of a Mott gap in 
the region U = 3 ... 4, which is caused by a change of the 
self-energy at the Dirac points of the Brillouin Zone from 
metallic to insulating behavior. This transition is in good 
agreement with finite-size extrapolations of variational 
QMC simulations and with continuous-time QMC calcu- 
lations based on cluster DMFT. As a result of short-range 
fluctuations, the critical Coulomb energy is significantly 
smaller than in single-site DMFT calculations. Also, a 
narrow pseudogap is found close to the Mott transition. 
Finally, the effective mass shows a moderate enhance- 
ment at finite U, but it does not diverge at the transition. 

The consistency between the ED and QMC calcula- 
tions for the honeycomb lattice, including the variation 
of the Mott gap with onsite Coulomb repulsion, suggests 
that, as long as the overall size of the Hilbert space is suf- 
ficiently large, yielding small enough level spacing, the 
use of only one bath level per impurity orbital can be 
adequate. This situation differs from the one for fewer 
sites or orbitals, where more bath levels per impurity 
level must be included to achieve sufficiently large Hilbert 
spaces. 
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